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Procedure  For  Identifying  Parameters 
Affecting  Storm  Runoff  Volumes 
In  A  Semiarid  Environment 


J.  R.  Simanton,  K.  G.  Renard,  and  N.  G.  Sutter  1 


Runoff  prediction  for  a  given  amount  of  precipitation 
has  been  the  primary  concern  of  many  hydrologic 
research  programs.  One  such  program  being  conducted 
by  the  Agricultural  Research  Service  on  the  Walnut 
Gulch  Experimental  Watershed  in  southeastern  Arizona 
has  resulted  in  several  different  mathematical  models  to 
predict  runoff.  These  runoff-predicting  models  deal 
mainly  with  precipitation  characteristics  and  the  effect 
of  these  characteristics  on  runoff  amount  and  peak  rate 
(3,  6).^  In  an  attempt  to  simplify  and  standardize  runoff 
prediction,  the  Soil  Conservation  Service  (SCS)  of  the 
United  States  Department  of  Agriculture  has  developed 
a  runoff-prediction  model.  The  SCS  National  Engineer- 
ing Handbook,  section  4,  Hydrology  (7)  (subsequently 
referred  to  simply  as  NEH-4),  illustrates  this  method  for 
predicting  watershed  runoff.  The  method  is  based  on 
known  precipitation  and  certain  watershed  character- 
istics, such  as  indices  of  the  soil-cover  complex  and 
antecedent  moisture  condition. 

The  hydrologic  soil  groups  are  divided  into  four 
classes  based  on  infiltration  and  soil-water  movement 
characteristics.  The  four  classes,  Usted  A  (the  lowest 
runoff  potential)  through  D  (the  highest  runoff  poten- 
tial), contain  each  soil  series  found  in  the  United  States. 

Cover  groups  are  determined  by  land  use  and  treat- 
ment classes.  Land  use  varies  from  row  crops  to  dirt 
roads  with  most  of  the  divisions  associated  with  various 
agricultural  uses.  Treatment  classes  consist  of  straight 
row  cropping,  contouring,  and  terracing.  The  land  use 
and  treatment  of  an  area  are  further  divided  into  three 
hydrologic  conditions  (poor,  fair,  or  good).  (See 
table  1.) 

Antecedent  moisture,  defined  as  the  summation  of 
the  5-day  precipitation  before  the  runoff-producing 
storm,  is  called  antecedent-moisture  condition  (AMC)  I, 
II,  or  III.  The  conditions  are  grouped  in  the  Handbook 
(NEH-4)  with  different  ranges  for  the  dormant  season 
and  for  the  growing  season  (table  2).  Because  all  runoff 
in  the  study  area  used  subsequently  in  this  paper  occurs 
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TABLE  I.— Runoff  curve  numbers  for  hydrologic 
soil-cover  complexes  for  watershed  condition  II 
and  I  J  =  0.2S 


Hydrologic  soO  group 

Land  use  Treatment  Hydrologic 

or  cover    or  practice    condition       A       B       C  D 


Pasture  or 


inge            ---  Poor  68  79  86  89 

Do            ---  Fair  49  69  79  84 

Do            ---  Good  39  61  74  80 

Do  Contoured  Poor  47  67  81  88 

Do             do  Fair  25  59  75  83 

Do            do  Good  6  35  70  79 


=  Initial  abstraction. 
S  =  Potential  maximum  retention  (inches). 

Source:  Soil  Conservation  Service  (7),  table  9.1,  p.  9.2 


TABLE  2.— Rainfall  groups  for  estimating  antecedent 
conditions  by  major  seasons 


Condition 

Five-day  antecedent  rainfall,  in  inches 

Dormant  season 

Growing  season 

I 

n 
III 

Less  than  0.5 
0.5  to  1.1 
Over  1.1 

Less  than  1.4 
1.4  to  2.1 
Over  2.1 

Source: 

Soil  Conservation  Service  (7),  table  4.2,  page  4.12. 

during  the  growing  season,  the  same  condition  ranges 
were  adapted  herein  as  follows:  AMC  I  =  0.00  inch  to 
1.40  inches;  AMC  II  =  1.40  inches  to  2.10  inches;  and 
AMC  III  is  over  2.10  inches. 

Knowing  these  watershed  conditions,  hydrologic  soil- 
cover  complexes,  and  antecedent  moisture,  one  can, 
through  the  use  of  tables  and  graphs  provided  in  the  SCS 
Hydrology  Handbook  (J),  determine  the  appropriate 
runoff  curve  number.  With  this  curve  number,  runoff 
can  be  predicted  for  a  given  watershed  rainfall  by 
graphical  techniques  or  by  solving  two  equations. 

The  sparse  runoff  data  available  in  the  western  United 
States,  when  the  original  NEH-4  was  developed,  Hmited 
the  verification  of  curve  numbers  for  conditions  found 
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in  the  semiarid  Southwest.  Studies  on  the  Walnut  Gulch 
Experimental  Watershed  are  suited  to  such  verification. 
These  studies  have  included  various  land  use  and 
management  practices  on  small  plots  and  small  natural 


vratersheds.  Curve  numbers  for  these  plots  and  small 
watersheds  are  determined,  and  a  method  is  presented  to 
extend  these  plot  curve  numbers  to  small  natural 
watersheds. 


DESCRIPTION  OF  STUDY  AREAS 


Runoff  data  were  used  in  this  study  from  five  areas  or 
groups  of  6-  by  12-foot  runoff  plots  on  the  Walnut 
Gulch  watershed  (fig.  1).  Tliese  groups  are:  KendaU 
plots,  located  at  the  upper  end  of  the  watershed  on  a 
grass-cover  site;  Tu-8  plots,  located  near  the  watershed 
center  on  a  grass-dominated  site  with  limited  shrub 
cover;  Lucky  Hills  and  Tu-9  plots,  located  in  the  north 
central  part  of  the  watershed  on  a  shrub-cover  site;  and 
the  Lamb's  Draw  plots,  near  the  watershed  outlet  on  a 
dense  shrub-cover  site.  The  soils  are  generally  shallow, 
with  below-average  infiltration  and  slopes  from  3  to  15 
percent. 

The  10  Kendall  plots,  with  their  20-percent  vegetative 
cover  and  moderate  grazing  practice,  would  be  classified 
as  fair  rangeland.  The  soil  at  these  plots  is  a  Hathaway 
gravelly  loam  with  8-  to  15-percent  slopes  and  is  in  the  B 
hydrologic  soil  group.  The  Hathaway  series  consists  of 
deep,  dark-colored,  well-drained  gravelly  medium  and 
moderately  coarse-textured  soils  over  very  gravelly 
coarse-textured  materials  at  moderate  depths  (I).  Four 
of  the  Kendall  plots  are  grazed  each  year,  while  the 


remaining  six  are  protected  from  grazing.  However, 
because  of  good  management  practices,  there  is  very 
little  difference  in  vegetative  basal  cover  between  the 
grazed  and  ungrazed  plots  (fig.  2). 

The  10  Lucky  Hills  plots  have  not  been  grazed  for  the 
past  8  years,  but  the  vegetative  cover  averages  only  30  to 
40  percent  so  they  would  be  classified  as  poor  rangeland. 
The  soU  at  these  plots  is  a  Laveen  gravelly  loam  with  3- 
to  8-percent  slope  and  is  in  the  C  hydrologic  soil  group. 
The  Laveen  series  consists  of  deep,  well-drained, 
medium-textured  soils  formed  in  old  calcareous  alluvium 
from  several  sources  (1).  No  rangeland  treatment  other 
than  ehminafing  grazing  has  been  performed  on  these 
plots. 

The  12  Tu-9  plots,  located  approximately  800  feet 
from  the  Lucky  Hills  plots,  have  a  vegetative  cover  that 
averages  40  percent  and  are  classified  as  poor  rangeland. 
Some  of  the  plots  would  be  classed  as  contoured  poor 
rangeland  because  of  treatments  imposed  on  them 
consisting  of  pitting,  clearing,  and  combinations  of  the 
two.  Each  treatment  is  rephcated  three  times.  As  a 
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Figure  2.— Kendall  6-  by  12-foot  grass-covered  (grazing  excluded) 
runoff  plot. 

check,  three  plots  were  left  untreated.  Plots  of  this 
group  are  located  on  a  Rillito  gravelly  loam  with  3-  to 
8-percent  slope  and  are  in  the  B  hydrologic  group.  The 
Rillito  series  consists  of  deep,  well-drained,  medium  and 
moderately  coarse-textured  gravelly  soils  formed  in 
calcareous  old  alluvium  (1).  Livestock  grazing  has  been 
excluded  from  these  plots  (fig.  3). 

The  four  Tu-8  plots  have  not  been  grazed  for  the  past 
8  years.  Vegetative  cover  averages  40  percent,  and  the 
plots  are  classed  as  poor  rangeland.  The  soil  at  these 
plots,  as  at  the  Kendall  plots,  is  a  Hathaway  gravelly 

RUNOFF  CURVE  NUM 

Runoff  curve  numbers  for  each  of  the  plot  locations 
were  determined  by  three  different  methods.  Method  1 
(table  and  graph  method)  involves  using  the  tables  and 
graphs  given  in  NEH-4.  Method  II  (graphical  method) 
uses  a  minimum  deviation  eye-fit  of  a  curve  on  the  graph 
of  actual  precipitation  versus  actual  runoff.  Method  III 
(mathematical  method)  uses  an  objective  minimum  least 
square  deviation  determination  using  a  digital  computer 
and  the  individual  precipitation  and  runoff  events.  Other 
methods  such  as  the  Phi  index  have  been  used  by  SCS 
but  were  not  included  in  this  study. 

Table  and  Graph  Method 

SCS  is  continuing  to  improve  its  methods  of  deter- 
mining runoff.  One  recent  improvement  in  the  pro- 
cedure has  been  a  figure  developed  especially  for  the 
hydrologic  conditions  experienced  in  this  area,  because 
they  realized  that  the  curve  numbers  determined  from 


Figure  3.— Tu-9  6-  by  12-foot  shrub-covered  (grazing  excluded) 
runoff  plot. 

loam  with  3-  to  8-percent  slopes  (1)  and  is  in  the  B 
hydrologic  group. 

The  12  Lamb's  Draw  plots,  also  classified  as  poor 
rangeland,  have  been  divided  into  three  treatment 
groups.  Four  plots  are  located  in  a  contour  pitted  area, 
four  in  a  root-plowed  and  seeded  area,  and  the  remaining 
four  in  a  control  area.  Vegetative  cover  averages  45 
percent  at  these  12  plots,  and  the  soil  is  the  same  as  the 
Tu-9  plots,  a  Rillito  gravelly  loam  with  3-  to  8-percent 
slope  (1). 

ER  DETERMINATION 

tables  1  and  2  produced  runoff  values  significantly  lower 
than  what  were  actually  observed  in  the  Basin  and  Range 
area  of  the  Western  United  States.  Figure  4,  which  was 
developed  from  figures  9.5  and  9.6  of  Chapter  9  in  the 
NEH-4  (7),  expresses  the  runoff  curve  number  as  a 
function  of  the  cover  density  and  hydrologic  soO  type 
for  various  vegetation  types. ^  The  procedure  for  deter- 
mining the  runoff  curve  number  from  figure  4  involves 
soil-cover  complex  associations.  For  the  Lucky  Hills  plot 
group,  the  vegetative  cover  averages  35  percent,  and  the 
soil  is  in  the  C  hydrologic  group.  The  curve  number  from 
the  Desert  Brush  curves  of  figure  3  associated  with  tliis 
soil-cover  complex  is  87.  This  procedure  was  followed 
for  the  remaining  plot  groups  with  the  exception  of  the 
Kendall  plot  group.  The  curve  number  for  this  group  was 
determined  from  the  Herbaceous  curves  of  figure  4. 
Table  3,  under  the  heading  "Method  I,"  lists  curve 
numbers  for  each  plot  group  determined  as  described. 

^Personal  communications  with  SCS. 
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Tlie  graphical  method  of  determining  runoff  curve 
numbers  was  from  the  actual  rainfall-runoff  data  ob- 
tained from  the  plots.  Rainfall  for  runoff  events  was 
determined  from  recording  rain  gages  located  adjacent  to 
each  plot  group.  Runoff  measured  in  the  55-gallon 
containers  was  converted  to  inches  runoff  from  the  plot. 
An  eye-fitted  curve  was  drawn  througli  the  plotted 
points  on  a  graph  of  rainfall  versus  runoff  for  each  plot 
group  and  treatment.  This  curve  was  then  compared 
with  the  curves  found  on  figure  5  and  then  assigned  the 
curve  number  from  the  figure  that  best  described  the 
developed  curve. 

An  example  of  this  method  for  the  data  from  Lucky 
Hills  plots  is  shown  in  figure  6.  Although  a  considerable 
scatter  is  experienced  in  the  measured  data,  changing 
curve  numbers  from  90  to  92  still  seems  to  fit  the  data 
fairly  well.  Therefore,  a  value  of  91  was  chosen  for  these 
data.  Included  in  figure  6  are  AMC  I  and  III  curves 
corresponding  to  the  average  antecedent  moisture  curve 
of  91 .  These  AMC  I  and  III  curves  were  determined  from 
table  10.1  of  NEH-4  (7).  These  two  curves  encompass 
nearly  every  observed  point  on  the  graph.  Column 
"Method  11"  of  table  3  lists  the  graphically  developed 
curve  numbers  for  each  plot  group  and  treatment. 

Mathematical  Method 

The  mathematical  method  of  determining  runoff 
curve  numbers  also  used  actual  rainfall-runoff  data.  The 
data  were  analyzed  using  a  digital  computer  program. 


TABLE  3.— Runoff  curve  numbers  developed  by  three  different  methods  for  the  Walnut  Gulch  Plots 

[Average  antecedent  moisture.] 

Plot  group  and  Number  o J  plots  Land  use  or  Method  I  Method  II  Method  III 

watershed  used  treatment 


Kendall  

6 

Fenced 

Do  . ,  

4 

Grazed 

Tu-8   

4 

Fenced  control 

10 

Fenced 

Tu-9   

3 

Fenced  control 

Do  

3 

Fenced  pitted 

Do  

3 

Fenced  cleared 

Do  

3 

Fenced  pit,  clear 

4 

Grazed  control 

Do  

4 

Grazed  pitted 

Do  

4 

Grazed,  root- 

plowed,  seed 

Kendall  Weir  #1  (143  acres) 

do 

Lucky  HiUs  Weir  M  (11  acres) 

Lucky  Hills  Weir  4^  (0.56  acre) 

do 

Watershed  63.004  (560  acres) 

79 

93 

93 

79 

94 

94 

81 

88 

90 

87 

91 

91 

81 

86 

85 

None  established 

87 

88 

do 

92 

92 

do 

89 

90 

81 

90 

86 

None  established 

82 

83 

None  established 

87 

90 

86 
87 
85 
85 
81 


Number  of  plots  used  to  determine  average  runoff. 


4 


RAINFALL  (P)   IN  INCHES 

Figure  5. -Runoff  prediction  from  curve  numbers  and  rainfall  amount. 


Given  a  set  of  actual  P  and  Q  values,  this  program  used 
an  iterative  "hill  climbing"  technique  to  solve  for  the 
optimum  value  of  S  in  the  SCS  runoff  equation: 

(p-a2s£ 

^     P  +  0.8S  ^  ' 

where  Q  =  Runoff  (inches) 
P  =  Rainfall  (inches) 

S  =  Potential  maximum  retention  (inches). 

This  equation  considers  the  initial  abstraction  to  be 
0.2S.  It  is  also  conceivable  that  additional  variance  of 
the  observed  versus  predicted  runoff  values  might  be 
explained  by  changing  the  initial  abstraction  term  in 
equation  1  from  the  0.2S  value  which  has  been  adapted 
by  SCS.  Although  such  an  optimization  was  not  in- 
cluded in  this  method,  it  could  be  easily  accomplished 
but  would  not  allow  direct  comparison  with  the  other 
information  in  NEH-4. 

Assuming  the  potential  maximum  retention  (5)  is  1 
inch  of  water  per  foot  of  gravelly  loam  soU,  the  initial 


abstraction  would  be  0.2  inch  if  only  1  foot  of  soil  is 
wetted,  which  is  a  phenomenon  generally  encountered 
with  the  short  duration  storms  such  as  are  considered 
herein.  Schreiber  and  Kincaid  (6)  have  shown  the  initial 
abstraction  to  be  0.25  inch  forThe  TU-8  and  TU-9  plot 
groups  which  closely  agrees  with  the  above  assumption. 
Given  a  set  of  P  and  Q  values,  the  least  sum  of  squares 
method  was  used  to  solve  for  some  value  f(S)  where 

n 

f(S)=  2   (QeJ-(Qoi)^and  (2) 
i  =  1 

Qej  is  the  estimated  runoff  using  a  trial  S  value  and 

an  observed  Pj 
Qoj  is  the  observed  runoff  corresponding  to  the 

observed  Pj  used  to  calculate  Qgj 
n     is  the  number  of  observations. 

The  initial  S  value  selected  was  3.00  inches  which 
corresponds  to  a  curve  number  of  around  77.  The  initial 
increment  of  potential  maximum  retention,  AS,  was 
then  set  to  be  0.5  of  the  initial  S.  f(S),  a  number,  is  then 
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Figure  6. -Lucky  Hills  plots,  precipitation  vs.  runoff,  1968-70  data. 


compared  with  f(S-AS).  If  f(S)  is  found  to  be  tlie 
smaller  of  the  two,  then  S  is  increased  by  AS.  and  the 
process  is  repeated.  Likewise,  if  f(S-AS)  is  smaller,  S  is 
decreased  by  AS,  and  the  process  repeated.  When  the 
condition 


f(S+AS)<f(S)<f(S+2AS) 


(3) 


is  met,  AS  is  negated  and  divided  by  two  since  the 
optimum  point  has  been  bypassed,  and  the  iteration 
process  is  repeated  in  the  opposite  direction  using  the 
smaller  increment  value  (fig.  7). 

This  technique  allows  "homing  in"  on  the  optimum  S 
value  by  bouncing  back  and  forth  across  the  boundary 
using  smaller  and  smaller  increments  for  AS.  When  AS 
becomes  less  than  0.005  inch,  the  process  is  terminated. 
Tliis  yields  an  approximate  accuracy  of  ±0.01  inch  for 
the  calculated  S  value. 


The  curve  number  for  any  set  of  P  and  Q  data  can  be 
determined  by  taking  the  calculated  S  value  and  apply- 
ing it  to  the  formula: 


CN  = 


1000 
10+S 


(4) 


Curve  numbers  obtained  for  the  plot  groups  and  their 
treatments  are  listed  in  Column  "Method  III"  of  table  3. 
The  calculated  curve  numbers  using  this  method  are  very 
similar  to  those  obtained  using  "Method  II."  The 
similarity  in  curve  numbers  developed  by  Methods  II  and 
III  indicates  that  either  of  these  methods  can  be  used 
successfully  to  determine  runoff  curve  numbers  from 
actual  rainfall-runoff  data.  For  many  field  locations 
without  a  digital  computer.  Method  II  is  the  most 
practical.  Method  III,  after  the  necessary  data  have  been 
digitized  on  computer  cards,  is  easier,  less  time  con- 
suming, more  objective,  and  hopefully  the  most 
accurate. 
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S  +  2AS 


PARAMETERS  AFFECTING  CURVE  NUMBERS 


Rangeland  Treatment 

Rangeland  treatment  effect  on  runoff  can  be  demon- 
strated by  differences  in  runoff  curve  numbers.  For 
example,  within  the  Lamb's  Draw  plot  group,  the  curve 
numbers  (Method  III)  range  from  83  for  the  pitted  plots 
to  90  for  the  root-plowed,  seeded  plots.  The  control 
plots  have  a  curve  number  of  86,  which  is  midway 
between  the  curve  numbers  developed  for  the  pitted  and 
root-plowed,  seeded  plots.  The  pitted  plots,  with  their 
greater  surface  water  storage  capacity,  would  be  ex- 
pected to  have  a  curve  number  lower  than  the  control 
plots.  The  curve  number  of  the  pitted  plots  should,  with 
time,  increase  as  the  pits  fill  mih  sediment.  The 
root-plowed,  seeded  plots,  with  their  initial  lack  of 
vegetative  cover,  might  be  expected  to  have  a  curve 
number  greater  than  the  control  plots.  The  curve 
number  for  the  root-plowed,  seeded  plots  should  de- 
crease as  the  seeded  grasses  improve  the  cover.  Ulti- 
mately, both  the  pitted  treatment  and  root-plowed  and 
seeded  treatment  may  be  expected  to  stabilize  at  lower 
curve  numbers  than  those  for  the  control  plots.  The 
Tu-9  plot  group,  with  its  treatments  of  pitting,  clearing, 
and  the  combination  of  the  two,  seems  to  be  nonre- 
ceptive  to  rangeland  treatments.  Curve  numbers  de- 
veloped for  each  treatment  were  higher  than  the  curve 
number  developed  for  the  control  plots.  Information 
from  the  Kendall  plot  group  suggests  that  proper  grazing 
intensity  may  not  alter  runoff  curve  numbers. 


Antecedent  Moisture  Conditions 

Although  the  frequency  distribution  of  antecedent 
moisture  was  not  uniform  within  the  moisture  classes,  an 
effort  was  made  to  determine  the  effect  of  antecedent 
moisture  on  the  runoff  curve  numbers.  An  AMC  of  less 
than  1.4  inches  existed  for  all  of  the  runoff-producing 
storms  analyzed  for  the  Walnut  Gulch  plots.  To  deter- 
mine the  effect  of  antecedent  moisture  on  curve 
numbers,  the  SCS  lowest  antecedent  moisture  clas-.  was 
divided  into  four  subconditions  defined  as  follows: 

O.OCKAMC  Ii<0.20<AMC  I2<0.50<AMC 


I3<0.90<AMC  l4<1.40 


(5) 


The  rainfall-runoff  data  for  each  plot  group  and  treat- 
ment were  then  grouped  into  these  four  classes  and 
curve  numbers  developed  using  Method  III.  Figure  7 
shows  the  change  in  curve  numbers  mth  each  change  in 
antecedent  moisture  subconditions.  The  curve  numbers 
for  the  Kendall  plots  changed  greatly  with  a  change  in 
antecedent  moisture.  The  curve  number  change  for  these 
plots  was  from  91  for  an  AMC  1 1  to  97  for  AMC  I4.  To 
demonstrate  the  importance  of  these  curve  number 
changes,  the  expected  runoff  from  an  assumed  1-inch 
rain  is  also  shown  in  fig.  8.  The  same  runoff  increase 
with  increased  antecedent  moisture  was  obtained  from 
the  Tu-9  and  Lucky  Hills  plot  groups  (fig.  8).  However, 
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Figure  8. -Influence  of  antecedent  moisi 

the  increase  was  not  as  large  as  it  was  for  the  Kendall 
plots.  The  decrease  in  curve  number  with  a  correspond- 
ing increase  in  AMC  for  the  Tu-8  plots  may  be  explained 
by  four  high -in  tensity  storms  occurring  in  the  AMC  1 2 
group.  Runoff  curve  number  changes  with  changes  in 
rainfall  intensity  are  discussed  later.  Limited  data  pre- 
vented developing  a  curve  number  for  AMC  I4  at  the 
Lucky  Hills  and  Tu-9  plots.  No  AMC  grouping  of  the 
Lamb's  Draw  rainfall-runoff  data  was  made  because  of 
the  short  period  of  record. 

Rainfall  Intensity 

For  the  semiarid  Southwest  and  its  characteristic 
precipitation,  the  maximum  15-minute  rainfall  depth  has 
been  shown  to  be  the  dominant  variable  in  determining 
runoff  from  small  areas  (3).  Three  intensity  classes  were 
selected  to  include  the  range  of  observed  maximum 
15-minute  rainfall  intensities  determined  for  each  storm. 
The  intensity  classes,  in  inches  per  hour,  were  defined 
as: 

0.25<Pi<1.00<P2<2.00<P3<4.40  (6) 


0.80r 


u. 


0.00 1-  1  1  ' 

AMC  AMC  AMC  AMC 

I  2  3  4 

ANTECEDENT   MOISTURE  GROUP 

—  0.20" 

—  0.50" 
0.90" 
1.40" 

ure  on  curve  numbers  and  observed  runoff. 

Method  III  was  again  used  to  determine  runoff  curve 
numbers  for  the  intensity  and  plot-grouped  rainfall- 
runoff  data. 

Curve  numbers  varied  appreciably  when  rainfall  in- 
tensity was  included.  An  increase  in  curve  number  for 
each  plot  group  occurred  with  each  increase  in  intensity 
class.  The  largest  curve  number  change  occurred  on  the 
Kendall  fenced  plots.  Here  the  curve  number  increased 
from  67  for  the  lowest  intensity  class,  Pij ,  to  92  for  the 
medium  intensity  class  Pj^  (fig.  9).  Relating  these  curve 
numbers  in  terms  of  runoff  from  a  1-inch  rainfall 
shows  that  an  intensity  class  of  (CN  =  67)  would 
produce  no  runoff,  whereas  an  intensity  class  of  Pj^ 
(CN  =  92)  would  produce  0.44  inch  of  runoff.  Limited 
data  for  the  Kendall  plot  group  prevented  developing  a 
curve  number  for  the  Pj^  intensity  class. 

The  Tu-8,  Tu-9,  and  Lucky  Hills  plot  groups  also 
showed  an  increase  in  curve  number  with  an  increase  in 
intensity  class  (fig.  9).  Rainfall-runoff  data  from  the 
Lamb's  Draw  plot  group  were  not  subclassed  because  of 
insufficient  field  data  for  the  2-year  record. 
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Figure  9. -Influence  of  rainfall  intensity  on  curve  numbers  and  observed  runoff. 


Combination  of  Antecedent  Moisture  and 
Rainfall  Intensity 

Rainfall-runoff  data  for  the  Lucky  Hills  plot  group 
were  then  grouped  into  both  antecedent  moisture  and 
intensity  classes.  This  plot  group  was  chosen  because  it 
has  the  greatest  number  of  rainfall-runoff  records.  Curve 
numbers  were  determined  for  each  combination  of 
intensity  and  antecedent  moisture  as  foUows: 

Plj-AMC  Ii,  Pij-AMC  I2,  Pij-AMC  I3,  Pj^-AMC  Ii, 

Plj-AMC  l2,Pi2-AMC  l2,Pi3-AMC  Ij. 

These  groupings  were  selected  because  there  was  a 
sufficient  number  of  storms  for  the  computations. 
Figure  10  shows  the  influence  of  AMC  and  rainfall 
intensity  on  the  curve  numbers  determined  by  Method 
III. 


The  curve  number  increases  shown  in  figure  9 
indicate  that  antecedent  moisture  in  some  soils  may  have 
a  greater  effect  on  runoff  at  low  rainfall  intensities  than 
at  high  rainfall  intensities.  This  may  explain  the  high 
curve  number  for  the  low  AMC  at  the  Tu-8  plot  group. 
This  plot  group  had  rainfall  intensities  of  greater  than 
2.00  inches  per  hour  for  50  percent  of  their  low  AMC 
storms.  Additional  data  need  to  be  gathered  before  this 
relationship  can  be  accurately  reported. 


Increasing  Drainage  Area 

The  precipitation  variabiHty  associated  with  runoff- 
producing  storms  in  the  semiarid  Southwest  is  a  limiting 
factor  in  curve  number  development  for  larger  water- 
sheds. To  correlate  rainfall  and  runoff  on  finite-size 
semiarid  watersheds,  Osborn,  Lane,  and  Hundley  (4) 
showed  that  three  evenly  spaced  rain  gages  are  needed 
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Figure  lO.-Influence  of  antecedent  moisture  and  rainfall 
intensity  on  curve  numbers  for  the  Lucky  Hills  plot 
group. 

for  a  1 -square-mile  watershed  and  that  a  network  of  five 
evenly  spaced  rain  gages  is  needed  for  a  10-square-mile 
watershed.  For  watersheds  larger  than  10  square  mOes, 
direct  correlation  between  rain  gage  networks  and  runoff 
decreases  rapidly.  Because  of  this  precipitation  vari- 
ability, a  distributed  model  with  flood  routing  should  be 
used  for  areas  larger  than  560  acres.  However,  as  has 
been  reported  for  larger  watersheds  in  the  semiarid 
Southwest,  surface  water  yield  per  unit-area  decreases 
with  increasing  drainage  area  (2,  5).  This  means  that  as 
the  watersheds  increase  in  size,  there  is  a  corresponding 
decrease  in  runoff  curve  number  due  mainly  to  the 
complex  drainage  systems  encountered  on  the  larger 
watersheds.  These  ephemeral  stream  channels  with  their 
coarse  alluvium  act  as  a  large  water  reservoir  and  account 
for  large  runoff  losses. 

Curve  numbers  for  rainfall-runoff  data  from  water- 
sheds larger  than  the  6-  by  12-foot  plots  were  developed 


by  using  Method  III.  Three  of  the  five  watersheds 
studied  are  located  near  the  Lucky  Hills  plot  group.  The 
soils  for  these  watersheds  are  classified  as  a  Rillito- 
Laveen  gravelly  loam,  and  the  vegetation  is  predomi- 
nantly shrubs  with  small  areas  of  grass.  The  watershed 
areas  are  0.56  acre,  11.00  acres,  and  560  acres.  The 
remaining  two  watersheds  are  located  near  the  Kendall 
plot  group.  The  soils  at  these  watersheds  are  classified  as 
a  Bernardino-Hathaway  gravelly  loam  with  a  vegetative 
cover  of  predominantly  grass.  The  areas  of  the  two 
watersheds  are  4.6  acres  and  143  acres. 

Curve  numbers  for  these  watersheds  are  listed  in  table 
3  under  Method  III.  As  would  be  expected,  the  curve 
numbers  for  these  larger  areas  are  lower  than  the  curve 
numbers  for  the  plot  groups.  From  these  curve  numbers, 
equations  were  developed  to  express  the  relationship 
between  watershed  size  and  curve  number.  Two  sets  of 
equations  were  developed  using  two  different  coeffi- 
cients (k)  in  the  equation: 


CN  =  kxe 


(7) 


where  CN  =  curve  number  for  drainage  area  x 
k  =  constant 
e  =  constant 

X  =  drainage  area  in  acres  for  which  CN  is 
being  determined. 

Three  equations  are  developed  in  each  set  to  show 
vegetation  difference.  One  equation  for  a  shrub-covered 
watershed,  one  for  a  grass-covered  watershed,  and  one 
for  the  combination  of  the  grass-  and  shrub-covered 
watersheds.  In  an  attempt  to  relate  a  plot  curve  number 
to  a  large  areas'  curve  number,  the  plot  curve  number 
was  used  as  the  k  value  in  the  first  set  of  equations  listed 
in  table  4.  As  is  shown  by  the  standard  error  of  the 
estimate  associated  with  the  equations,  this  method  is 
not  very  accurate  in  predicting  the  larger  area's  curve 
number.  The  k  value  for  the  second  set  of  equations  is 
derived  from  the  least  sum  of  squares  of  the  hne  drawn 
through  the  plotted  points  on  a  graph  of  curve  number 
versus  drainage  area.  This  second  set  of  equations  and 
their  standard  error  of  estimate  are  also  listed  in  table  4. 
The  standard  error  of  estimate  for  this  set  of  equations  is 
considerably  less  than  the  standard  error  of  estimate  for 
the  first  set  of  equations.  These  equations  can  be  used  to 
determine  curve  numbers  for  areas  smaUer  than  560 
acres  with  similar  soil-cover  complexes  as  those  for 
which  the  equations  were  developed. 
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TABLE  4.— Prediction  equations  showing  the  effect  of  drainage  area  on  curve  numbers 


Set  of  equations  and  type  watershed  Prediction  equation  Standard  error  of  curve  number 


Shrub-covered   CN  =  91x~0-0088  +10.6 

Grass-covered   CN  =  93x"0-0078  4g  j 

Combined  grass  and  shrub   CN  =  92x"0-°08°  ±14.3 

Second  set: 

Shrub-covered   CN  =  85.75  x""-""°Z  ±1.6 

Grass-covered   CN  =  88.00x""-°^°^  ±1.1 

Combined  grass  and  shrub   CN  =  86.74x~0-00*^8  ±3.5 


\\Tiere  x  =  drainage  area  in  acres  with  the  largest  watershed  tested  =  560  aaes. 


CONCLUSIONS 


Soil  Conservation  Service  design  runoff  curve  num- 
bers developed  from  actual  rainfall-runoff  data  for  the 
Walnut  Gulch  runoff  plots  are  somewhat  higher  than 
those  recommended  by  the  Soil  Conservation  Service's 
Hydrology  Handbook.  This  is  partly  due  to  the  lack  of 
channels  on  the  plots  which  on  larger  watersheds  cause  a 
decrease  in  runoff  amount  because  of  channel  abstrac- 
tions. Cur\'e  numbers  developed  for  areas  larger  than  the 
6-  by  12-foot  runoff  plots,  but  with  very  similar 
soil-cover  complexes,  decrease  with  increasing  drainage 
area.  Equations  developed  to  express  this  relationship 
have  been  presented  but  are  limited  to  watersheds  with  a 
similar  soil-cover  complex  and  less  than  560  acres.  The 
higher  than  recommended  curve  numbers  developed  for 
the  plots  are  also  due  to  the  high  rainfall  intensity 
associated  with  runoff-producing  storms  in  the  semiarid 
Southwest.  This  high  rainfall  intensity  is  a  dominant 
factor  in  runoff  production  from  semiarid  watersheds 
and  should  be  considered,  along  with  soil-cover  complex 
and  antecedent  moisture  when  developing  a  runoff  curve 


number  for  a  watershed  in  cUmatic  provinces  like 
southeastern  Arizona.  The  antecedent  moisture  condi- 
tion classes  as  listed  in  NEH-4  are  insufficient  to  describe 
the  antecedent  moisture  conditions  ranges  commonly 
found  in  the  semiarid  Southwest.  A  finer  breakdown  of 
these  antecedent  moisture  classes  is  necessary  to  show 
curve  number  changes  due  to  antecedent  moisture. 

The  development  of  runoff  curve  numbers  from 
rainfall-runoff  data  for  watersheds  larger  than  1  square 
mile  is  very  difficult  because  of  nonuniform  precipita- 
tion distribution. 

For  accurate  runoff  prediction  on  small  watersheds  in 
the  semiarid  Southwest,  the  SCS  runoff  curve  number 
method  would  better  include  the  factors  of  rainfall 
intensity,  Hmited  areal  extent  of  the  rainfall,  and 
watershed  drainage  area,  in  addition  to  the  soil  and 
vegetation  parameters.  Runoff  design  from  larger  water- 
sheds could  then  be  included  by  summing  the  individual 
subwatersheds  runoff  estimates  and  reducing  the  com- 
bined estimate  for  transmission  losses  in  major  channels. 
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